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FOREWORD 


This  report  describes  the  sonar  signal  characteristics  and  classification 
research  study  performed  by  Honeywell  Inc.  ,  Systems  and  Research  Center, 
under  U.  S,  Naval  Ship  Systems  Command  Contract  No.  N00024-68-C- 1219. 
Mr.  G.  Miller  of  the  Naval  Ship  Systems  Command  provided  direction  for 
the  program. 

Principal  authors  are  Dr.  Duane  H.  Tack  and  Dr.  Eugene  E.  Yore,  of 
Honeywell's  S&RC  Research  Department.  The  authors  gatefully  acknowledge 
the  technical  contributions  and  leadership  of  Richard  M.  Powell,  principal 
investigator,  and  later  project  supervisor,  Charles  Johnson  was  program 
manager. 


The  Hone3’well  number  assigned  to  this  report  is  12120-FR. 
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ABSTRACT 


(C)  A  study  of  the  space -time  properties  of  active  sonar  echoes  is  made  with 
a  view  toward  target  classification  on  the  basis  of  shape  and  aspect  at  long 
range  in  the  bottom  bounce  mode.  Spectral  and  time  cross-correlations  of 
the  echoes  at  two  space  points  are  examined  for  feasibility  of  signal  proces¬ 
sing  with  respect  to  requirements  on  the  transmitted  signal  and  transmitter- 
target -receiver  geometries.  STARLITE  processing  (spectral  correlation) 
is  the  most  promising  of  the  two,  but  a  detailed  examination  of  the  echo 
properties  leads  to  more  severe  restrictions  on  the  transmitted  signal  and 
receiver  separations  than  previously  indicated.  Ideally,  a  wide,  flat,  spec¬ 
tral  window  (created  by  the  transmitted  signal),  within  which  a  band-limited 
spectral  correlation  is  performed,  is  desired.  Filtering  and  spreading  of 
the  signal  spectrum  upon  reflection  off  an  irregular,  stratified,  impedance 
bottom  distorts  the  spectral  window.  This  places  additional  constraints  on 
the  operational  parameters  for  which  the  space-time  properties  of  the  echo 
yield  target  information  beyond  that  which  is  available  at  a  single  receiver. 
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SECTION  I 
INTRODUCTION 


(C)  Classification  clues  obtainable  from  time-processed  echoes  that  are  use 
ful  in  short-range  sonar  operation  become  somewhat  less  auspicious  at  long 
range  in  the  bottom  bounce  mode  [3].  Accurate  extraction  of  echo  length  is 
more. difficult  in  the  low  signal-to-noise  ratio  environment  created  by  long 
propagation  paths  and  bottom  bounce  losses.  Target  localization  within  the 
beam  makes  target  shape  and  aspect  information  unattainable  from  time- 
processed  echoes.  As  a  result,  the  other  clues  such  as  positional  consis¬ 
tency,  localization  within  the  beam,  and  Doppler  become  more  significant  in 
long-range,  time-processing  sonars  [3],  However,  space-time  processing 
renews  the  possibility  of  extracting  target  shape  and  aspect  [l,  2],  Since 
this  information  is  a  valuable  addition  to  the  other  classification  clues  avail¬ 
able  at  long  range,  this  .report  examines  the  space-time  properties  of  sound 
fields  reflected  from  targets  characterized  as  linear  arrays  of  point  reflec¬ 
tors  with  a  view  toward  extracting  target  shape  and  aspect  information  for 
classification  at  long  range  in  the  bottom  bounce  mode, 

(C)  Space-time  analysis  of  the  target- reflected  sound  field  suggests  either 
spectral  or  time  cross-correlations  for  classification  signal  processing. 
Spectral  cross-correlations  have  been  considered  previou.sly  in  STARLITE 
[2],  However,  the  received  signals  at  two  space  points  do  not  precisely 
possess  the  properties  attributed  to  them  in  the  STARLITE  analysis.  As  a 
result,  a  reexamination  of  the  STARLITE  concept  was  required,  leading  to  a 
special  type  of  spectral  cross -correlation  and  somewhat  more  severe  con¬ 
straints  on  its  applicability.  For  time  cross-correlations,  the  space-lime 
properties  of  the  target- reflected  sound  field  require  a  compression  of  the 
time  scale  of  one  of  the  received  signals  at  two  space  points.  Time  cross - 
correlations  for  target  classification  therefore  require  more  complicated 
signal  processing  than  spectral  correlations. 
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ASSUMPTIONS 

(C)  The  basic  purpose  of  this  study  was  to  examine  the  effects  of  bottom 
reflections  upon  the  received  sound  field  and  their  ramifications  to  signal 
processing  for  classification.  In  order  to  limit  the  scope  of  such  a  task  to  a 
manageable  level,  extraction  of  target  shape  and  aspect  in  the  bottom  bounce 
mode  has  been  chosen  for  investigation,  partly  for  the  reasons  given  above 
and  partly  because  it  is  one  of  the  more  interesting  long-range  classification 
clues  from  the  standpoint  of  the  space-time  correlation  properties  of  the 
received  bottom  bounce  echoes.  To  simplify  the  problem  and  concentrate  on 
the  salient  aspects  of  bottom  bounce  effects  on  target  shape  and  aspect  classi¬ 
fication,  stationary  targets  and  transmitter- receivers  are  assumed,  trans¬ 
missions  through  water  are  taken  to  be  ideal,  and  ambient  noise  and  rever¬ 
beration  are  neglected.  Also,  multipaths  involving  surface  reflections  are 
signored,  and  as  a  result,  target  and  transmitter- receiver  depths  are  not  a 
factor  in  the  analysis.  Ocean  depth  enters  only  through  the  range  of  incident 
angles  for  bottom  reflections  since  loss  of  signal  power  over  the  propagation 
path  is  immaterial  in  the  assumed  noiseless  medium. 

(U)  The  ocean  bottom  is  best  characterized  by  an  irregular,  stratified  impe¬ 
dance  boundary.  To  model  these  characteristics  as  manifested  in  the  bottom 
reflection  coefficient,  a  physical  approach  using  the  Kirchoff  radiation 
formula  is  taken  to  model  bottom  reflections  [7].  Because  reflection  and 
first-order  scattering  are  linear  processes,  the  effects  of  bottom  reflections 
can  ultimately  be  modeled  by  a  sum  of  random,  attenuat  ag,  delay  line  filters. 
The  spectral  properties  of  these  filters  depend  upon  the  angle  of  incidence  of 
the  signal,  the  amount  of  irregularity  in  the  sediment  layers,  and  the  values 
of  the  acoustic  parameters  (density,  velocity,  attenuation  rate)  characteri¬ 
zing  the  sediment  layers. 

(U)  Time  and  scope  limitations  on  this  contract  did  not  permit  a  detailed 
investigation  of  the  effect  of  bottom  reflections  on  target  classification  on 
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(U)  the  basis  of  shape  and  aspect.  However,  preliminary  estimates  based 
on  the  spectral  properties  of  the  bottom  transfer  function  indicate  that  bottom 
bounce  STARLITE  may  be  practical  at  long  range  where  incident  angles  are 
near  grazing  and  the  effects  of  surface  roughness  and  penetration  of  the  boun¬ 
dary  by  the  sound  wave  minimized. 


SUMMARY 

(C)  Section  II  examines  the  general  properties  of  the  target -refle^’ted  sound 
field  and  the  resulting  implications  to  space-time  processing  for  classifica¬ 
tion  by  spectral  and  time  cross-correlations.  The  analysis  leads  to  some¬ 
what  more  severe  geom.etric  and  signal  bandwidth  constraints  for  spectral 
processing  than  indicated  by  STARLITE  [2j,  However,  the  limitations  of 
time  cross-correlations  are  even  more  severe;  hence,  STARLITE -like 
signal  processing  appears  to  be  the  most  useful  for  shape  and  aspect 
classification  by  space-time  analysis. 

(U)  Section  HI  considers  the  effects  of  bottom  reflections  and  scattering  in 
the  specular  direction  upon  the  spectral  crsoss -correlations.  The  model  for 
the  bottom  reflection  coefficient  is  presented  in  Appendix  A. 


(U)  Sectio."'  IV  covers  conclusions  and  recommendations. 
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SECTION  U 

SPACE-TIME  ANALYSIS 


(C)  Thii.  sertion  examines  the  general  space -time  properties  of  the  sound 
field  1‘efiected  from  a  target  array  of  point  scatterers  and  the  resulting  impli¬ 
cations  to  signal  processing  for  target  classification  on  the  basis  of  shape  and 
aspect.  Two  types  of  signal  processing  are  considered:  cross-correlation  of 
the  spectra  of  signals  received  at  two  space  points,  i.e. ,  ST  Ail  LITE -like 
processing  [2],  and  time  cross-correlations.  It  will  be  seen  that  the  received 
signals  do  not  precisely  possess  all  the  properties  attributed  to  them  in  the 
STARLITE  analysis,  with  the  result  that  somewhat  more  severe  restrictions 
are  placed  on  the  transmitted  signal  spectrum  and  transmitter-target-receiver 
geometries  for  successful  extraction  of  target  shape  and  aspect. 


GENERAL  SPACE-TIME  PROPERTIES  OF  REFLECTED  SOUND  FIELD 


(U)  For  the  purposes  of  deriving  the  spectral  properties  of  the  target-reflected 
sound  field,  assume  that  an  arbitrary  array  of  (M  +  1)  point  reflectors  of 
scattering  amplitude  a^  (i  =  0,  1,  . . . ,  M)  is  insonified  by  an  omnidirectional 
transmitter  through  a  nonrefracting,  nonreflecting,  dispersionless  medium. 

In  the  far  field,  the  received  signal  spectrum  at  the  j  ‘  (omnidirectional)  re¬ 
ceiver  is  then  simply  , 

(i)  ,  (i) ' 


M 

w'  (u)  =  s^{u)  r. 

J  i=o 


a.e 

1 


J 


(1) 


where  is  the  sound  velocity,  Srp(u)  is  the  transmitted  signal  spectrum, 

"  is  the  distance  from  the  transmitter  and  r.^^^  is  the  distance  from  the 


T 


'} 
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receiver,  respectively,  to  the  i  “  point  reflector  in  the  target.  At  long  range 


^  r,j,(o)  + 


where  L-  is  the  distance  between  the  0*^  and  i*^  reflectors,  L  =  0,  rrT,^°^ 

,  \  l  /■»  O  /-x  1 

and  r^  '  are  the  nominal  path  lengths  to  the  target,  and  '  and  a}  ^  are 
the  direction  cosines  of  reflector  i  relative  to  the  specular  directions  of  the 
transmitter  and  receiver,  respecti\’-ely,  in  coordinate  systems  having  their 
origins  at  the  0**^  reflector.  Since  r^°^  »  at  long  range,  target  shape 

and  aspect  information  is  essentially  contained  only  in  the  phase  of  the  re¬ 
ceived  signals.  Thus,  approximating  r^^^  by  r^^^  in  the  amplitude  but  retain¬ 
ing  Equation  (2)  in  the  phase,  the  spectra  of  the  received  signals  at  two  differ¬ 
ent  space  points  (j  =  1,  2)  can  be  written  in  the  form 


Wj(w) 


where 


W^(u))  =  H2(u)S^(w)e'j‘^S 


TT  /  \  Mu 

H.(u)  =  s  b.e  '' 
i=o 


iVl 

H2(w)  =  b'  2  b^e 


are  the  target  transfer  functions  for  the  two  viewing  angles,  and  where 
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b" 

t. 

J 


H 


1 


ri^^VrgCo) 

r  +  r 
T  ^^3 


T  1 


(Xr 


(i) 


Tir 


O', 


+  a 


+  a 


J 


(5) 


Although  Equations  (3),  (4),  and  (5)  have  been  derived  under  idealized  condi¬ 
tions,  the  basic  functional  forms  of  W  ^  and  are  unchanged  in  a  refracting, 
reflecting,  snd  dispersive  medium  so  long  as  refraction  and  dispersion  are 
negligible  over  the  length  of  the  target  and  there  are  no  reflecting  surfaces 
within  distances  of  the  order  of  the  transmitted  signal  duration  times,  c^, 
from  the  target  in  the  specular  directions  for  transmission  and  reception. 
Under  these  conditions,  refraction  affects  only  the  nominal  path  delays,  t., 

J 

multipath  returns  can  be  resolved  in  time,  and  dispersion  and  reflection  alter 

only  the  shape  of  the  transmitted  signal  spectrum,  S,j,(u)  (requiring  that  Srp  be 

replaced  by  S.(io)  in  W  c.f.  Section  III),  leaving  the  target  functions,  H.(u), 

J  3  3 

unchanged.  Thus  more  generally  we  can  write 


Wj(u)  =  H^(u)Sj(u)e'j“*l  =  Wj(u)e"j‘*^*l  ( 
W2'(w)  =  H2(u)S2(u)e~j'^^2  =  W2(w)e'j‘^S 


(6) 
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(C)  In  the  case  of  linear  arrays  typical  of  submarine  targets,  =  a„, 

j  J 


=  a.,  and  hence,  =  e  for  all  i.  Under  this  condition 


H2(w)  =  b'Hj(eu) 


(7) 


and  it  is  recognized  that  H2((<))  represents  a  compression  (e  >  1)  or  expansion 
(e<  1)  of  the  frequency  scale  of  Hj(«)  upon  change  in  aspect  angle.  It  is  this 
property  of  the  target  transfer  functions  as  viewed  from  two  different  space 
points  which  is  exploitable  in  extraction  of  target  shape  and  aspect.  The  func¬ 
tion  of  signal  processing  for  classification  is  then  to  determine  whether  or  not 

Equation  (7)  holds,  and  if  so,  estimate  e  to  determine  target  aspect  angles 
T 

Of  and  and  finally,  estimate  Tj^j  to  determine  target  length  Lj^.  However, 
signal  processing  for  classification  must  be  accomplished  with  the  received 
signal  spectra,  and  it  is  evident  that 

W^(u))  ^  h'  W^(eu)  (8) 

even  when  the  path  delays,  t^,  are  removed  by  appropriate  shifts  of  the  time 
scales  of  the  two  received  signals.  Basically,  this  results  from  the  fact  that 
the  transmitted  signal  does  not  undergo  a  corresponding  time  scale  expansion 
or  compression  by  the  factor,  e,  upon  reflection  off  an  individual  scatterer  on 
the  target.  Thus,  the  inequality  of  (8)  must  be  overcome  by  appropriate  de¬ 
sign  of  the  transmitted  signal  and  the  signal  processing  methods  in  order  that 
the  equality  of  (7)  can  be  observed,  when  true. 


SPECTRAL  CROSS-CORRELATIONS 

(C)  STARLITE  [2]  was  the  first  space-time  classification  scheme  to  use  the 
properties  of  the  transfer  functions  of  linear  targets  expressed  in  Equation  (7 ). 
In  particular,  it  was  recognized  in  STARLITE  that  the  inherent  periodicity  in 
the  transfer  function  of  a  line  target  varies  with  aspect  angle,  and,  over  a 
narrow  frequency  band,  tie  difference  in  period  between  the  target  transfer 
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(C)  functions  viewed  at  two  slightly  different  aspect  angles  can  be  approxi¬ 
mated  by  a  relative  shift  in  the  frequency  spectra  of  and  Hg.  As  a  result, 
the  STARLITE  classification  scheme  is  to: 


1)  Compare  the  frequency  spectra  of  two  return  signals  to  see  if 
they  are  the  same  except  for  a  shift  in  frequency.  If  so,  the 
target  is  a  line  target  and  one  can  proceed  to  step  2. 


2)  Measure  the  frequency  shift  to  estimate  aspect  angle. 

3)  Measure  the  least  period  of  H.  to  estimate  target  length, 

3 


All  three  of  the  above  tasks  can  be  accomplished  by  a  cross-correlation  of 
the  frequency  spectra  of  the  two  received  signals. 


(C)  Two  major  constraints  on  the  geometries  for  which  STARLITE  is  applic¬ 
able  as  listed  in  reference  [2]  are  the  so-called  correlation  condition 


Af 


B 


w 


<0.5 


(9) 


where  Af  is  the  (cyclic)  frequency  shift  and  is  the  bandwidth  of  the  trans¬ 
mitted  signals,  and  the  resolution  condition 


>  1 


(10) 


These  conditions  were  arrived  at  through  consideration  of  the  properties  of 
and  H2.  but  apparently  ignoring  the  realities  of  inequality  (8).  Therefore, 
it  is  necessary  to  reexamine  the  STARLITE  conditions  and  the  special  form 
of  the  spectral  cross-correlation  required  in  view  of  this  inequality. 
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Spectral  Correlation  Function 


(C)  Since  and  Wg  do  not  possess  the  properties  of  and  H2  for  all  fre¬ 
quencies,  it  is  necessary  to  concentrate  on  portions  of  their  frequency  spectra 
that  do,  i.e. ,  narrow  bands  where  the  are  constant,  or  approximately 
so.  As  pointed  out  in  the  STARLITE  analysis  [2],  the  transmitted  signal 
essentially  provides  a  spectral  window  within  which  it  is  possible  to  view  the 
target  properties  contained  in  and  Hg*  However,  in  attempting  to  measure 
these  properties  of  the  target  from  and  W2>  the  shape  of  the  spectral  win¬ 
dow  must  not  be  permitted  to  obscure  or  influence  them.  To  this  end,  we 

consider  the  spectral  cross-correlation  function 

.  A 


r  ° 

\ 

1  A 

W 

'^2 

1  ' 

u  -  0 

1 

/ 

(11) 


where  is  the  (radian)  carrier  frequency  of  the  transmitted  signal,  Q  is  the 
relative  (radian)  frequency  shift  introduced  in  the  two  spectra,  the  asterisk 
denotes  the  complex  conjugate,  and  Re  indicates  that  the  real  part  of  the  inte¬ 
gral  is  taken. 


(C)  Ideally,  the  signal  spectra  should  be  constant  (S^  =  A^)  in  the  integration 
band  [w^-  -  A/2,  +  A/2],  for  then  Wj(w)  =  Hj(w)  in  the  integral.  If  the  band 

A  is  small  and  the  change  in  target  aspect  angle  between  the  two  receivers  is 
small,  so  that  1 1  -  c  |  «  1,  for  the  integration  interval  we  can  approximate 


=  2  ttN.  +^.  +  (u  -  u^)t^ 

*  2nN.  +  y.  +  (u  -  u^)eTj^  (12) 

-  2nN.  +  Xi  +  (w  -  u^)ti 
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(C)  where  N.  is  an  integer.  Under  these  conditions 


K,2(q) 


M  M  sin^(t.  -tj  r 

=  A  A  b'AE  sb.b  -  2 - icos  §(Ti+  Tj)  "  (Ci '  Xj> 

i=0  3*o  J  ^ 


Spectral  Correlation  Constraints 


(C)  The  diagonal  terms,  j  =  i»  in  Equation  (13)  have  maxima  at 


(Ci  -  Xi) 


^  ^  =  (1  -  e)  (0^;  i  =  1.  2. . . . .  M 


by  subtraction  of  the  second  line  of  Equation  (12)  from  the  first  line  and 
evaluating  the  difference  at  u  =  This  is  the  desired  frequency  shift  for 
line  targets  that  we  seek.  They  also  have  maxima  at  n  "  Xi^  *  Znin, 

m  =  1,  2.  ....  leading  to  an  ambiguity  with  which  we  will  deal  later.  The  off- 
diagonal  terms  do  not  have  maxima  at  Q.  hence  they  must  be  rendered  negli¬ 
gible  [in  order  that  the  maximum  of  o::cur  at  o  =  Q]  by  imposing  the 

condition 


I"  ■'min'  '  V 


if  they  are  to  contribute  less  than  10  percent  of  corresponding  diagonal  com- 

A 

ponents  to  the  value  of  R 

(C)  Consider  next  the  effect  of  nonflat  signal  spectra  on  R^2‘  illus¬ 

trative  simplicity,  assume  that  the  band-pass  signal  spectra  can  be  approxi¬ 
mated  by  S.(u)  A. cos  P(u  -  u^)  in  the  neighborhood  of  the  carrier  frequency. 
Then  the  effect  of  such  on  R^2  replace  the  sin  x/x  term  in  Equation (13) 

by  the  term 


C. .=  T  cos  /3n 
ij  ^ 


sin|(T.  -T.)  ^  <Ti  -  Tj-23)^  sin  -|(t.  -  Tj  +  23) 

I  (xi  -  V  Lf  f  (n-j"23) 
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(C)  For  C. .  to  be  insensitive  to  n,  ci^.«  1  is  required.  If  we  take  u  ±#  to 
*3  •  o  ^ 

coincide  with  the  half -power  points  of  S^(u),  |3  *  n/2A  and  differs  by  less 
than  5  percent  of  its  value  at  Q  =  0  if 


_Q 

A 


S  0.  25 


(16) 


The  values  of 


C.. 

13 


,  i  7^  j,  are  again  rendered  negligible  by  condition  (15). 


(C)  It  is  necessary  that  have  a  single  global  maximum  within  the  range 

of  frequency  shifts  to  be  encountered  if  a  unique  q  is  to  be  obtained  from  the 
correlation.  A  unique  global  maximum  can  be  ensured  by  requiring  that 


IOtJ  s  2n  (17) 

since  0<Tj^<...<  and  there  can  then  be  only  one  q  for  which  every 
diagonal  element  in  Rj^2(n)  Is  maximum.  Furthermore,  Q  is  then  readily  iden¬ 
tifiable  as  the  first  global  maximum  encountered  either  to  the  left  or  right  of 
Q  *  0.  As  a  counter  example,  consider  the  case  of  just  three  scatterers 
(M  =  2)  and  T2  "  3x^/2.  If  condition  (17)  were  violated  and  the  desired  global 
maximum  occurred  at  0  =  2n+  0  <  y  s:  2tt  .  there  is  also  a  global  maxi¬ 
mum  atQTj  =  y"2n=Ci"\j^‘'^TT>  since  then  0  72  ®  C2  "  X2  "  ^ 

is  closer  to  the  origin  than  Q  »  it  would  be  mistaken  for  the  desired  frequency 
shift  and  an  erroneous  estimate  of  g,  and  hence  target  aspect  and  length, 
would  be  made.  Although  this  example  is  based  upon  commensurate  target 
highlight  delays  in  practice  they  would  probably  only  have  to  be  approxi¬ 
mately  commensurate  to  generate  this  kind  of  error  if  the  ambiguity  condi¬ 
tion  (17)is  not  met,  since  the  maxima  of  the  fluctuating  components  in 
are  quite  broad.  Note  that  conditions  (17)  and  (15)  yield  a  somewhat  stronger 
constraint  than  (16),  since  together  they  imply  12tt/A  s  ^  |nTj|52rr, 

or  I^yAl  1/6.  Thus,  (15)  is  a  stronger  condition  for  correlation  than  (16). 
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(C)  Limitations  on  the  minimum  measuratile  frequency  shift  in 
more  a  function  of  the  ability  to  determine  that  a  global  maximum  of  R^^g 
been  achieved  in  the  presence  of  noise  than  to  measure  the  period  of  the  most 
rapidly  fluctuating  component  in  implied  by  the  STARLITE  resolution 

condition  (10).  At  signal-to-noise  ratios  sufficient  for  detection,  it  is  reason¬ 
able  to  expect  that  one  would  be  able  to  ascertain  the  difference  between  a  zero 
and  a  maximum  in  each  significant  correlation  component  in  R^g*  Since  the 
term  with  the  smallest  delay  requires  the  largest  variation  in  Q  between  a 
zero  and  an  adjacent  maximum,  the  resolution  condition  for  determination  of 
a  global  maximum  of  Rj^g  presence  of  noise  is  approximately 


QtiI  s  2 


(18) 


(C)  In  summary,  the  significant  constraints  on  spectral  processing  for  classi¬ 
fication  of  arrays  of  point  scatterers  are: 


1)  Correlation  condition 
"rminl  ^ 


2)  Ambiguity  condition 
'l 


In  T,  Is  2tt 


3)  Resolution  condition 


n  Ti 


TT 


Additional  STARLITE  conditions  [2]  are  the  Fresnel  condition,  which  is  not 
constraining  at  long  range,  and  the  thickness  condition  which  limits  the  ap¬ 
parent  length-to-diameter,  D^,  of  submarine  targets  by  the  condition  [2] 


sin  0 

- 15 

s 


>4 
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(C)  where  <ji  is  the  target  aspect  angle  (see  Figure  1).  These  constraints  can 
be  untenable  if  significant  scatterers  in  the  target  are  too  close  together.  For 
example,  the  correlation  condition  is  violated  for  target  highlight  delays  sepa¬ 
rated  by  less  than  10  msec  at  a  signal  bandwidth,  B  ,  of  600  Hz.  It  would 
seem,  therefore,  that  successful  classification  of  targets. that  behave  acous¬ 
tically  as  arrays  of  point  scatterers  by  spectral  processing  cannot  depend  on 
extraction  of  all  target  highlight  delays,  but  rather  must  rely  on  extraction  of 
only  those  widely  spaced  highlights  such  as  the  bow,  conning  tower,  and  stern 
of  a. submarine. 


Example 


(C)  To  examine  briefly  the  implications  of  the  constraints  of  conditions  (15), 
(17),  and  (18)  on  reciever  separations  and  the  transmitted  signal  spectrum, 
consider  the  simple  two-dimensional  STARLITE  geometry  of  Figure  1  where 
the  transmitter  and  one  receiver  are  positioned  side  by  side.  For  this  geom¬ 
etry,  the  aspect  angle,  0,  is  defined  to  lie  in  the  range  of  180®,  with 

beam  aspect  at  0°  or  180°  and  bow  or  stern  aspect  at  90° ,  Also,  the  following 
approximations 


€  - 


D 


[1  +  cos  Y  -  cot  0  sin  Y  ] 
(O)^ 


_  Y 


cot  0 


(19) 


sin  0 


(o) 


<  <1, 


(o). 


=  rj^(o)  rg 


(o) 


Y  <<  1).  Assuming 


hold  at  long  range  (D/r,j, 
a  target  length  of  90  m  and  just  three  significant  scattering  surfaces,  "  12 
(0  =  tt/2)  =  120  msec  for  c^  =  1.5  Km/ sec.  It  is  also  reasonable  to  assume 
that  T2y3  ^  ^''"2/3  target  delays  associated  with  the  conning  tower. 

Signal  bandwidth  requirements  can  be  fixed  using  the  STARLITE  thickness 
constraint  [2]  which  limits  the  minimum  aspect  angle  for  classification  to 
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Figure  1.  STAR  LITE  Geometry 
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(C)  0  =-  15°  in  this  case.  Then,  since  12^3 

satisfied  for  bandwidths  B  =  A/2n.  i.e., 

w 


^  T 

nain 


^  ^2/2' 


constraint  (15)  is 


400  Hz  5  B  s  600  Hz 
w 


(20) 


Assuming  that  signal  bandwidths  of  one-tenththe  carrier  frequency  can  be  readily 
generated,  let  f^  =  =  6  KHz.  Maximum  receiver  separations,  D,  are 

required  when  is  a  minimum.  Thus,  for  =  '’*2/3'  constraints  (17)  and 
(18)  evaluated  at  q  =  (1  -  e)  yield 


2.08  X  10“^  s  f  Icos  0l  s  8.33  x  10“^ 


(21) 


under  these  conditions.  Acceptable  receiver  separations  at  ranges  r  = 
of  20K  and  40K  yards  and  at  6  =  90°  implied  by  (21)  are  plotted  in  Figure  2  as 
a  function  of  aspect  angle,  0.  Corresponding  STAJRLITE  constraints  [2],  c.f. 
Equations  (9)  and  (10),  are  about  4/3  larger  on  the  minimum  separation,  and 
for  all  practical  purposes,  unbounded  above  for  the  conditions  of  this  example. 
A  receiver  separation  of  about  175  yards  would  permit  extraction  of  target 
shape  and  aspect  for  15°  s  0  ^  60°  and  120°  0  s  165°  at  0  =  90°  and  20K  ^  r  g 

40K  yards.  Although  this  is  probably  an  acceptable  receiver  separation  dis¬ 
tance,  the  constraints  are  tight  and  do  not  permit  much  flexibility  in  range 
and  bearing  for  fixed  D.  In  addition,  the  signal  bandwidth  requirepients  are 
an  order  of  magnitude  larger  than  the  bandwidths  used  in  the  initial  STAR  LITE 
trials  [2].  It  therefore  appears  that  some  flexibility  in  choosing  receiver 
separations  as  function  of  range  and  bearing  is  necessary  for  target  shape  and 
aspect  classification  by  spectrail  cross-correlation. 


TIME  CROSS-CORRELATIONS 

(U)  After  compensation  for  the  path  delays,  t.,  the  functional  form  of  the 

1 

received  signals  in  the  time  domain  is 
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TARGET  ASPECT  ANGlE,  .  ,  DEGREES 


(C)  P^igure  2,  Receiver  Separation  Constraints  versus 
Target  Aspect  Angle  for  the  Example  (U) 
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Wi(t)  *  2  b.3^(t  -  n)  (22) 

M 

Wgit)  =  b'  2  b.s  (t  -  e  tJ 

In  order  that  and  Wg  achieve  maximum  correlation  when  e-  =  &  a  com¬ 
pression  (e  >  1)  or  expansion  (s  <  1)  of  the  time  scale  by  the  amount  e  must 
be  performed  in  order  that  the  signal  component  reflected  off  the  i^^  scatterer 
in  have  the  same  delay  as  the  corresponding  component  in  Wj  for  every  i. 
Thus,  we  seek  the  time  scaling  e  of  w-  such  that 

A 

M 

W2{t)  -  b'  2^  b^S2[e  (t  -  T^)]  (23) 

Suppose  that  an  estimate  of  e  is  g,  yielding 

M 

W2(t)=b'  2  b^S2[e  (t  -  uT.)]  (24) 

i=o  ^ 

where  ^  =  e  /c  ,  and  consider  the  time  cross-correlation  of  w.  and  w 

1  2 


Ri2('i)  =J  T)dt  =  2^:^  J  Wj(u))W2  dw 


Since 


~  b'  /  i  M 

W„(u)  =  —  2  b.e'J'^^'^i 


~  2  ~  “"r 

e  '  e  '  1=0 


R^2^t)  can  be  written  in  the  form 
M  „ 

BjjW  -b'  b,  B^^3^[T-U-„)Ti] 
M  i- 1  i 


M  1=  i 


where 
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1  " 

Kg  g  (t)  =  —  r  (28) 

1  2  2rrG  J  *  ?  ' 

“OD 

(U)  The  functional  form  of  Rj2tT)»  Equation  (27),  is  sufficient  to  give  us  some 
clues  as  to  the  utility  of  the  time  cross-correlation  in  space -time  classifica¬ 
tion  on  the  basis  of  target  shape  and  aspect.  Although  the  carrier  frequency 
has  not  been  explicitly  removed  in  the  preceding  analysis,  we  are  interested 
in  the  correlation  of  the  envelopes  of  and  W2,  or  alternatively,  Ihe  envelope 
of  R ^^2^1  )• 

(U)  At  e  *  ei  u  ■  1  and  the  classification  condition  is  achieved  by  maximizing 
the  diagonal  elements  in  Rj^2^t)  at  t  =  0.  However,  the  off-diagonal  elements 
do  not  have  maxima  at  j  =  0.  Hence,  they  must  be  rendered  negligible  by  a 
signal  bandwidth  constraint  similar  to  that  of  Equation  (14). 


Resolution  Condition 

(C)  Resolution  of  the  diagonal  components  in  places  even  more  severe 

restrictions  on  its  properties.  It  we  require  a  variation  of  the  slowest  vary¬ 
ing  component  in  the  envelope  of  R^2(t)  of  the  same  order  of  magnitude  as 
that  which  was  required  for  the  slowest  varying  component  in  Rj^2^0)»  then  we 
require  that  R  [(1  -  g)  j.]  be  small  compared  with  R  (0),  or  approxi- 


I  1  -  e  I  ^  1  (29) 

Comparison  of  this  resolution  condition  with  that  of  (18)  for  spectral  correla- 

A 

tions  at  n  =  n  -  (1  ■  e  )<Jq  shows  that  spectral  processing  holds  a  distinct  ad¬ 
vantage  over  time  processing  in  that  the  carrier  frequency,  f^,  rather  than 
the  signal  bandwidth,  B  ,  appears  on  the  left  in  the  resolution  condition  for 
spectral  processing.  Since  |1  -  s  |  is  directly  proportional  to  the  receiver 
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(C)  separation  and  >  B^,  the  resolution  condition  for  spectral  correlations 
will  be  met  with  shorter  receiver  separations  than  those  required  for  time 
correlations.  As  an  example,  for  the  same  conditions  as  the  example  dis¬ 
cussed  in  connection  with  spectral  cross-correlations,  condition  (29)  reduces 
to 

I Y  cos  0  I  >  0.  08  (30) 

which,  by  comparison  with  condition  (21),  indicates  that  more  than  an  order 
of  magnitude  larger  receiver  separations,  D,  would  be  required  for  time 
correlations  than  for  spectral  correlations.  The  receiver  separation  could 
be  reduced  by  increased  signzil  bandwidths,  but  it  is  doubtful  that  an  order  of 
magnitude  reduction  could  be  achieved  and  still  keep  the  carrier  frequency  in 
the  low-to-mid-kilocycle  range  required  for  long-range  propagation. 


Linear  FM  Signals 


(U)  Before  concluding  the  discussion  of  time  correlation,  the  important 
special  case  of  linear  FM  (LFMj signals  should  be  examined,  since  it  turns 
out  that  the  time  cross-correlations  of  the  envelopes  of  compressed  LFM 
pulses  have  the  resolution  properties  associated  here  with  time  processing, 
whereas  direct  time  cross- correlations  of  the  envelopes  of  the  received  LFM 
echoes  have  the  properties  of  spectral  correlations. 


(C)  Consider  received  LFM  pulse  components  s.  in  w.  of  the  form 

3  3 

3(Wf^t  +  —  pt^) 

s,(t)=A,  e  ^  ^  [u(t+  Tq)  -  u(t  -  Tp)] 


3 


3 


(31) 


where  u(t)  is  the  unit  step  function.  For  correlation  after  pulse  compression 
we  assume  that  the  received  signals  are  processed  by  a  matched  filter  having 
the  transfer  function 


Hj(w)  =  exp 


j<“  -  “o* 

2p 


2  -\ 


(32) 
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(C)  with  the  result  that  the  compressed  signal  pulse  is  [22] 


s/(t)  =  A, 
J  3 


/2|iT, 


sin  ^Tjjt 

"W 


j(Wnt 


1  9  TT 

+  T) 


From  Equation  (28),  the  time  correlation  of  the  envelopes  of  s'j  and  s^ 
(at  e  =  1)  is  then 


^s'  s'  “  ■^1^2  — 5 
®  1®  2  ^  ^  2tt^ 


sin  H'TqT 


But  =  M-Tq/tt  for  the  LFM  pulse;  hence,  for  Rg/  g/ 


to  have  its  first  cor- 


W  •  U  ‘  S  iBo 

relation  zero  at  T  =  tt/h^Tq  <  |1  -  e  |tj^.  condition  (29)  is  required.  We  note, 
therefore,  that  reduction  of  the  reverberation  noise  level  in  STAR  LITE  pro¬ 
cessing  by  LFM  pulse  compression  as  suggested  by  Wiekhorst  [23]  is  ob¬ 
tained  only  at  the  expense  of  increased  receiver  separation  requirements. 


(C)  It  is  well  known  that  the  envelope  of  an  LFM  pulse  is  a  good  approximation 
to  its  spectrum;  hence,  a  time  cross-correlation  of  the  envelopes  of  uncom¬ 
pressed  LFM  echoes  should  be  a  good  approximation  to  their  spectral  cross¬ 
correlation.  To  investigate  this  case  mathematically,  consider  that  within 
the  time  interval  [Tq  +  T^  -  ['Tq  +  ‘  whichever 

is  the  smaller,  the  received  signals  can  be  written  in  the  form 


Wj(t)  =  s  b.s^(t  -  Tj^)  =  s^(t)  £ 

1  S 


W2(t)  =  b'  £  b.S2(t  -  ei:)  =  S2(t)b  '  £  b.  e 


u  e-3H^  ^^t)T.  -|(iTi^] 


=  A^  s^(t)Zj(t) 
-j[(WQ  +  lAt)€T.  - 
^Ag  ^S2(t)  Z2(t) 


when  Sj  is  given  by  Equation  (31),  and  is  therefore  the  envelope  of  w^  within 
the  specified  time  interval.  To  the  same  order  of  approximation  as  in  the 
spectral  correlations  [c.f.  Equation  (12)]  we  can  write 
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(Wq  Ti  =  2nN.  +  Ci 
(Wq  +  Jit)  ct^  =“  2  rrN.  +  X.  + 


(36) 


Then,  a  time  correlation  analogous  to  the  spectral  correlation  of  Equation  (11) 
is  T 

R^2<cr)  =“  Re  r  Zfi-  (t  +  J)  dt 

-T 


M  M 


sin  P-T(t-  -  T.) 

.iL 


aTA^A^b'E^  E^bjb,^T(Ti-Tj) 

-  -  V  -  ^ 


COS 


Jia 

Y(Ti+Ti)  (37) 


where  T  ^  -  t  ~  a  is  chosen  to  ensure  integration  over  a  time  inter- 

0  m  "^max  ^ 

val  where  all  (M  +  1)  echo  components  in  w^  and  vf^  are  non-zero. 


(C)  Identifying  Q=  ji-a  and  A=  2^iT  =-  2  ttB^,  if  >>T^  +  ^max’  ^12^^^ 
exactly  the  form  of  Equation  (13),  except  for  an  additional  term  in 

the  argument  of  the  cosine  function.  Thus,  in  addition  to  the  constraints 
outlined  for  R^^(fi)  in  the  discussion  of  spectral  cross-correlations,  there  is 
the  condition  that 


IL 

2 


-  Xi  |,  all  i 


(38) 


to  be  imposed  if  used  as  an  approximation  to  the  spectral  correla¬ 
tion,  Dividing  Equation  (38)  by  using  Equation  (14),  tJ-  =  ttB^/Tq, 

and  approximating  1 1  +  e  ]  2,  this  condition  reduces  to 


<< 


(39) 


for  the  most  stringent  case,  i  =  M.  Condition  (39)  indicates  that  the  ratio  of 
signal  bandwidth  to  carrier  frequency  should  be  much  less  than  the  ratio  of 
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(C)  signal  duration  to  maximum  target  highlight  delay,  a  condition  which  is 
easy  to  meet  with  LFM  pulses.  Condition  (39)  is  offered  as  an  alternative  to 
the  so-called  "transient"  condition  of  [23]  which  seems  to  be  based  on  re- 
quiring  that  (  1/2  |  s  2tt  ,  or  B^/(2Tq)|  s  1,  in  this  notation. 

Though  condition  (39)  is  less  restrictive,  it  is  felt  to  be  a  more  realistic 
indication  of  the  errors  involved  in  approximating  the  spectral  correlation, 
by  the  time  correlation,  Rj2(ff)* 

(C)  Because  larger  receiver  separations  and  more  complex  signal  processing 
(time  compression  of  one  received  signal)  are  required  for  time  correlations, 
it  is  felt  that  the  spectral  correlation  method  is  more  suitable  to  space-time 
processing  for  target  classifications  on  the  basis  of  shape  and  aspect  at  long- 
range  in  the  bottom  bounce  mode.  We  include  time  correlations  of  the  envel¬ 
opes  of  uncompressed  LFM  echoes  as  a  special  form  of  spectral  processing 
in  the  above  statement.  Therefore,  Section  III  of  the  report  is  devoted  to 
examining  the  effects  of  bottom  reflections  and  scattering  in  the  specular 
direction  upon  the  spectral  processing  for  classification. 
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SECTION  ni 

EFFECT  OF  BOTTOM  BOUNCE  ON  TARGET  CLASSIFICATION 


FORMULATION  OF  THE  BOTTOM  BOUNCE  TARGET 
CLASSIFICATION  PROBLEM 

(C)  As  the  discussion  of  Section  II  has  shown,  it  is  desirable  that  the  fre¬ 
quency  spectium  of  the  signal  be  flat  within  the  integration  interval  chosen 
for  the  frequency  cross- correlation  of  the  received  signal  spectra.  Pre¬ 
sumably,  the  transmitted  signal  spectrum,  S,p(u ).  will  be  chosen  to  have 
properties  which  approximate  this  condition.  However,  in  propatating  through 
a  random  inhomogeneous  medium  and  reflecting  off  rough  impedance  bound¬ 
aries,  the  signal  spectrum  undergoes  filtering  and  spreading  before  arriving 
at  the  receiver.  It  is  our  purpose  here  to  examine  the  nature  of  the  filtering 
and  spreading  that  the  signal  spectrum  experiences  upon  reflection  from  the 
bottom,  and  its  effect  upon  our  ability  to  extract  the  desired  target  charac¬ 
teristics  when  operating  in  the  bottom  bounce  mode. 

(C)  There  are  essentially  two  main  approaches  to  the  theory  of  reflection, 
reverberation,  and  scattering  in  random  inhomogeneous  media:  (1)  The 
physical  approach  as  discussed  by  Chernov  [4],  Tatarski  [5],  and  Laval, 
et  al.  [6],  in  relation  to  a  statistically  homogeneous  medium  characterized 
by  a  random  refractive  index,  and  by  Tolstoy  and  Clay  [7]  in  relation  to 
reflections  from  irregular  surfaces,  and  (2)  the  phenomenological  approach 
as  discussed  by  OHshevskii  [8]  and  Middleton  [9],  for  example,  wherein  the 
inhomogenieties  are  introduced  as  a  random  distribution  of  point  scatterers 
in  the  medium.  The  phenomenological  approach  has  the  advantage  of  sim¬ 
plicity  in  handling  complex  geometries  and  the  higher  order  statistics  re¬ 
quired  for  system  analysis  and  design.  However,  it  suffers  from  the  inability 
to  predict  the  impulse  response  (or  equivalently,  the  transfer  function)  of  the 


23 


CONFIDENTIAL 


CONFIDENTIAL 


(C)  collection  of  scatterers.  This  must  be  supplied  at  some  point  by  com¬ 
parison  of  the  theory  with  experiment.  But  the  success  or  failure  of  sonar 
signal  classification  in  the  bottom  bounce  mode  by  STARIiTE  techniques  de¬ 
pends  on  the  spectral  properties  of  the  bottom  transfer  functions  within  the 
integration  band  selected  for  the  cross- correlation  of  the  frequency  spectra 
of  two  spatially-separated  received  signals.  We  therefore  choose  to  take  the 
physical  approach,  since  the  basic  functional  form  of  the  bottom  transfer 
function  can  be  determined  from  measurable  physical  properties  of  bottom 
sediments  as  manifested  in  the  reflection  coefficient  [4]. 


(U)  To  formulate  the  bottom  bounce  problem,  we  follow  Tcistoy  and  Clay  [7] 
and  use  the  Kirchoff  radiation  formula  of  classical  diffraction  theory  which 
relates  the  pressure  interior  to  a  bounded  medium  to  the  pressure  on  the 
boundary  and  the  Green's  function.  Let  the  far-field  acoustic  pressure  at  a 


distance  r^  from  the  transmitter  be  given  by 
Pt  ^^T  ”  r 


where  is  the  source  strength,  d^  is  the  source  directivity,  =  w/c^  is 
the  wave  number  for  propagation  at  velocity  c^,  and  harmonic  time  dependence 
(~e^*^)  is  assumed.  Referring  to  Figure  3,  the  outgoing  pressure  wave  at  the 
point,  resulting  from  reflection  and  scattering  of  p.p  off  the  rough 

boundary,  z  =  0,  is  approximately  [10] 
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Figure  3. 


Geometry  for  Acoustic  Paths 
Bounce  Mode 


in  the  Bottom 
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where  =  Xrpi  +  y,pj  are  the  coordinates  of  the  bottom  plane  relative  to  the 
principal  axes  of  the  projection,  drpC'^rj.)>  of  the  source  directivity  function 
on  the  plane,  z  =  0,  and 


(i) 


YT 


(i) 


k  = 


Ic 

T 


h  ~  ^Tx  ^  ^Ty 


(42) 


is  a  vector  equal  to  the  difference  in  direction  cosines  of  the  vector  wave 
numbers  for  incident  and  reflected  waves,  is  the  random  eleva¬ 
tion  of  the  rough  bottom  relative  to  the  plane,  z  =  0,  and  w)  is  the 

bottom  reflection  coefficient  appropriate  to  the  illuminated  area,  drp(Xrj,). 


For  targets  located  far  from  the  reflection  point.  Equation  (41)  can  be  appr'^xi- 
mated  by  assuming  that  all  deviations  in  and  over  the  transmitter 
bottom  aperture,  d_(Xrr.),  from  their  values  in  the  specular  direction  are 
negligible,  i.e.,  0rp(^)  =  0  and  =  2  cos  Yp..  In  addition,  the  ' 
approximated  by 


can  be 


(i) 


(o) 


+  L.  oc,^ 


(i) 


(43) 


where  is  the  distance  between  the  first  (i  «  0)  and  (i  +  l)  (denoted  by  the 
superscript  i)  point  reflector  on  the  target,  =  0,  and  is  the  direction 
cosine  of  reflector  i  relative  to  the  specular  direction  of  the  bottom  reflected 
signal  in  a  coordinate  system  having  its  origin  at  the  first  point  reflector. 


Retaining  Equation  (43)  in  the  phase,  but  approximating  r^, 

(i) 

“Jtv  lirp  ’  Irp  T  i_,  M 

(i) 


(i) 


r. 


(o) 


in  the 


amplitude.  Equation  (41)  reduces  to  ,  , 

-jk  [r  +  r„'°'  +  L.  O' 

-jp^cosT,j,k^e 


p(r’n 


r  rp 


(O) 


/ 


— * 

T  X^  w)g 


■JkoYTC^ 


(44) 
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Recognizing  that  the  return  signal  emanating  from  the  i  point  reflector 
of  scattering  amplitude,  a^^,  at  a  distance,  r,  from  the  scatterer  is 


Pi(r)  = 


a  second  application  of  the  Kirchoff  radiation  formula  to  the  same  order  of 

th 

approximation  as  Equation  (44)  yields  for  the  pressure  at  the  j  "  receiver 


p/rj)  = 


- ^jdx^d.(x,)A.(?.,a,)e 


where  subscript  j  on  a  variable  denotes  the  same  type  of  function  as  defined 

til 

for  Equation  (44),  but  now  as  applies  to  the  j  receiver. 


Received  Signpl  Spectra 


Finally,  since  the  response  of  a  linear  system  to  a  simple  harmonic 

wave  is  the  transfer  function  of  the  system,  given  a  transmitted  signal,  S:j.(t), 

with  frequency  spectrum,  SrTp(u),  the  frequency  spectrum,  W.  tw),  of  the 

th  J 

received  signal  at  the  j  receiver  is  determined  by  Equations  (46)  and  (44), 
with  p  replaced  by  S~(ar).  Writing  W.  '(w)  so  as  to  explicitly  indicate  its 

>  0  i.  J 

frequency  dependence 


Wj'(u)  =  w^Sy(u) 


'  J:;£-,{dT(XT)d/VAT(Xx.“)VXj.w)e  (47) 
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(C)  we  can  write 


M  -UT. 
«!<“>  =  i%o 


M  T 

H„(w)  =  b'L  b.e 
-s  1=0  1 


by  defining 


b.  =  b 
1  1 


(1) 


/-  u  (2)  /.  (1)  ^2^ 

■  i  /  i  "  (0)  r~ 

r^cosij 


.  all  i 


T.  =  T . 
1  1 


(1) 


t/2)  + 

^  =  « V™liT  >  i  =  1'  2,  ....  M 

i  "t  "l 


Also,  the  signal  spectra,  S.(w),  can  be  written  in  the  form 

J 


S  (u)  =  G.(u)S^((o) 

J  J 


(51) 


(52) 


(53) 


where 


G.(u)  =  u  cC-  £■ 

J  Xnr 


J.  w)Aj(Xj^  w)e  ^  ^  ^  ^  |  (54) 


is  the  transfer  function  associated  with  the  bottom  reflections  encountered  by 

til 

S,p(w)  in  traveling  to  the  j  receiver.  It  is  seen  by  Equations  (50)  and  (51) 
that  the  signal  spectra,  Wj^(u)  and  W2(u),  have  the  desired  form  for  performing 
target  classification  ala  STA.RLITE,  provided  Sj^(w)  and  S2(u)  are  nearly 
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(C)  constant  over  the  integration  interval  of  Since  S^(b))  is  presuma¬ 

bly  chosen  to  be  sufficiently  flat  over  the  interval,  success  or  failure  of  the 


classification  scheme  ultimately  rests  on  the  properties  of  G.(u). 


Therefore, 


we  direct  our  attention  to  the  frequency  dependence  of  these  transfer  functions 


which  behave  like  random  filters  in  the  communication  channel. 


Expectation  of 


(C)  Because  of  the  roughness  of  the  bottom  surface  and  randomness  in  the 
reflection  coefficients,  the  statistical  expectation  of  is  of  interest. 

For  reasons  presented  later  in  this  discussion,  the  frequency  cross-correla¬ 


tion  for  bottom  bounce  operation  is  taken  in  the  form 

(Jq  +  2 

<Rj2(q)>  .  <Be  r  Wj  u.-  f)w;  y+f  d^> 

J  A  ' 


“o  ■  2 


A 

•"o  ■  2 


<G^(u  -  ^ 

where  the  bracket  <>  denotes  the  statistical  average  over  all  appropriate 
random  variables.  The  expectation  <GjG2  >  includes  both  c  he  rent  and  inco¬ 
herent  scattering  in  the  specular  direction  from  the  bottom  sediments.  No 
effort  is  made  to  separate  the  two  here,  since  it  is  the  properties  of  <R]^2^ 
which  are  of  interest,  not  those  of  G  ^  and  G2  in  and  of  themselves.  It  is  to 
be  noted  from  Equation  (55)  that  loss  of  coherence  in  the  specular  direction 
by  the  signal  upon  reflections  from  the  boundary  has  no  effect,  per  se,  on  the 
target  classification  scheme.  All  we  seek  is  a  proper  vehicle  for  creation  of 
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(C)  a  spectral  window  within  which  to  view  the  properties  of  -  n  /2) 

H,  (w  +  o/2).  Hence,  the  coherence  properties  of  the  reflected  signal,  S.(w), 
relative  to  the  transmitted  signal,  S,j,((o),  are  immaterial.  What  does  matter 
is  that  and  S2  have  sufficient  coherent  energy  within  the  band  [u^  -  A/2  , 

+  A/2]  to  raise  above  the  reverberant  and  ambient  noise  level  and 

that  the  distribution  of  that  energy  within  the  band  be  sufficiently  flat  so  as 
not  to  obscure  the  target  properties  contained  in  and  Hg. 


(U)  At  ranges  of  the  order  of  20  to  40  thousand  yards,  one  can  reasonably 
expect  essentially  complete  overlap  of  the  bottom  illumination  functions  d.^,, 
dj,  and  d2.  As  a  result,  fourth-order  statistical  moments  of  the  reflection 
coefficients  and  surface  roughness  are  required  in  computing  <  Gj^G2  >  .  This 
overlap  is  desirable,  since  a  high  degree  of  correlation  between  Sj^(w)  and 
S,,((jj)  improves  the  prospects  for  successful  target  classification.  However, 
this  presents  complications  relative  to  performing  the  spatial  integrations 
over  the  bottom  aperture  functions  when  weighted  by  these  higher-order  mo- 
ments.  To  reduce  the  complexity  of  <GjG2  >  ,  and  yet  preserve  the  essential 
sameness  of  and  S„,  it  will  be  assumed  that  there  is  no  statistical  depen- 
dence  of  the  reflection  coefficients  and  surface  roughness  between  the  out¬ 
going  and  return  reflections.  Physically,  this  corresponds  to  a  situation  where 
there  is  a  large  separation  between  the  transmitter  and  two  relatively  closely- 
spaced  receivers,  so  that  there  is  complete  overlap  of  d^  and  d2.  but  none 
between  d,j,and  either  d^  or  d2.  In  addition,  it  will  be  assumed  that  the  surface 
roughness  is  statistically  independent  of  the  reflection  coefficient.  The  expec- 
tation  <G^G2  >  is  then 


<Gj(u  -  ^)G2(u  +  ^)  >  *  (u  -  J  )£  -  At^T'  “  '  z' 


')  A  *  (y^',  u  +  2  )  > 


<exp  [-  ^  Cx^Xt*  '  2^  Cx<Xt'^^  ^  ^ 


.  ,  2  a  ,, 

(U  -  A  )* 


'R  '‘^r'|‘^R*Xr>'*h*Xr'*  ‘=^Ap(Yj,.u- 


iAp'(xR'-  “  +  2)> 


(56) 


<exp  {-  ■  2^  Cr^Xr)  ■  2  ^ 
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where  it  is  assumed  that  =  dg  =  dp,  =  r\2~  Tip  and  ^2  ~ 

Under  these  assumptions,  the  expectation  of  the  bottom  transfer  functions 
can  be  written  as 


<G^(u  -  §  )  G2*((J  +  (w«  n)  Ip(w, 


(57) 


where  I,p  and  Ip  have  the  same  functional  form  if  d^  and  dp,  A,p  and  Ap,  and 
rp  and  Tp  each  have  similar  functional  forms  and  distributions.  Assuming 
this  also  is  the  case,  the  problem  is  reduced  to  consideration  of  the  integral 


I(w,  Q )  =  I  -  T 


2> 


d(x)d(x')  <  A 
•  <expi-jri 


X-u  -2 


A* 


-*  /  , ,  j.  n 
X  .  w+  ^ 


(58) 


Evaluation  of  I(u,  n)  requires  specification  of  the  bottom  illumination  function, 
d(x),  the  functional  form  of  the  reflection  coefficient.  A,  and  the  joint  proba¬ 
bility  distributions  of  the  surface  roughness  and  random  parameters  in  a*  As 
described  in  Appendix  A,  the  bottom  reflection  coefficient,  A.  is  assumed  to 
have  the  form  ^ 

N 

A(X.“)  =  S  (59) 

n=o 


where  X.  is  a  reflection  coefficient  associated  with  the  interface  between  two 
n 

layers  in  the  stratified  medium,  is  an  inverse  velocity  associated  with  the 

wave  propagation  in  a  given  layer,  and  is  the  layer  thickness.  Based  on 

measurements  of  attenuation  in  bottom  sediments,  X.  and  0  are  taken  to  be 

n  n 

complex  constants.  The  randomness  in  A  is  assumed  to  be  characterized  by 
irregular  layer  thicknesses,  solely. 


The  details  of  the  evaluation  of  I(w,  n)  under  the  assumptions  of  a  gaussicin 
illumination  function  and  joint  gaussian  distributions  for  Q  and  i,  are  carried 
out  in  Appendix  B.  The  result  for  these  and  other  assumptions  made  in  the 
course  of  the  analysis  is  that 
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I(w.  n)  =“  -j  )e 


-  ^  (u^  +  J  ) 


where 


N  N 

L  L  X  J  (w,  n) 
n=o  m=o 


-  j(u  -I-  j  (G.  +1  )i*^»  -  i  (u  -  f  )2  K„  -  ‘  (u.  +  I 


'*■'  ■  ,?o 


n  n 

K„  =  L  L  (3  0  a 

q=o  p=o  P  ^P 


n)  =  (2TTDxDy)' 


2tt  DxDyLx  Ly 

nm  •'nm 


F  J  - 

Ei  u 


U  4 


In  addition,  Z  is  the  mean  of  -t,  ,  ct  is  the  covariance  of  I  (v)  and  l  (v) , 

2  9  _  q  qp  q  ^  'p'x' 

is  the  variance  of  Q{y),  Dx  and  Dy  are  characteristic  lengths  associated 

with  the  principal  axes  of  d(y),  Lx^^  and  Ly^^  are  correlation  lengths  asso¬ 
ciated  with  a  composite  spatial  correlation  coefficient,  xO,  defined  by 

Equations  (B15)  and  (B17),  along  the  principal  axes  of  drp(x),  is  a  coef¬ 
ficient  associated  with  the  transformation  of  the  correlation  lengths  appro¬ 
priate  to  the  principal  axes  of  ^  „  into  Lx  „  and  Ly  E.(x)  is  the  expo- 
nential  integral  as  defined  in  Equation  (B25),  is  Euler’s  constant,  and 
finally 


2  2  2 
.  »  r,0^ 


min(n,  m) 
L 

q,p=o 


a  ^  Re(|3  3  ) 

qp  q  p 
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SPECTRAL  PROPERTIES  OF  BOTTOM  TRANSFER  FUNCTION 

(C)  The  ultimate  purpose  of  this  investigation  vvould  be  to  examine 
in  the  vicinity  of  its  maximum  at  a  frequency  shift,  fi,  in  order  to  estimate 
how  much  deviation  from  the  desired  condition,  o  =  (l-e)u)^,  is  encountered 
in  the  bottom  bounce  mode  of  operation  when  the  target  is  a  line  target.  In 
addition,  one  would  like  to  design  the  transmitted  signal  spectrum,  S,p((o) ,  so 
as  to  compensate  for  errors  in  0  introduced  by  fcrttom  reflections.  Because 
of  the  complexity  in  the  integral  in  Equation  (55)  and  the  uncertainty  in  the 
acoustic  parameters  characterizing  the  bottom  sediments  over  a  variety  of 
bottom  types,  such  an  investigation  requires  an  extensive  computer  simulation 
of  the  bottom  bounce  classification  problem.  Unfortunately,  a  computer  sim¬ 
ulation  of  this  magnitude  was  not  within  the  scope  of  the  present  contract. 


Smooth  Boundaries 

(C)  Some  light  can  be  shed  on  the  possibilities  for  success  of  the  STAR  LITE 
classification  scheme  in  the  bottom  bounce  mode  by  examining  the  frequency 
dependence  of  I(tj,  Q)  in  the  vicinity  of  the  carrier  frequency,  since, 
ideally,  I(w,  n)  would  be  flat  over  the  interval  [w^-  A/2,  +  A/2]  and  inde¬ 

pendent  of  the  frequency  shift,  Q.  For  a  nearly  smooth,  regular,  stratified 

impedance  boundary,  such  that  ut  «  1,  coherent  reflections  dominate 

o  nm 

and 


1{U,  n)-(2nDxDy)7u  - 


where  the  real  and  imaginary  parts  of  the  exponent  have  been  exposed  by 
writing 
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(C)  Consider  first  the  diagonal  terms  in  the  double  sum  of  Equation  (64).  For 

1  p 

m  =  n,  the  exponent  in  the  addend  reduces  to  (-2uii.  +  jn|Ji  ).  Removal  of 

i  i  n  n 

the  term  1/2  n  ([a  -  K  ),  for  the  diagonal  terms  at  least,  is  the  reason  why 
n  m 

the  average  of  -  q/2)  W2  (w  +  n/2)  was  chosen  for  the  cross-correlation 
of  the  frequency  spectra,  since  this  term  potentially  generates  a  strong  ampli¬ 
tude  dependence  on  Q  if  an  unsymmetric  frequency  shift,  such  as  in  W^(u)  -Q) 

Wg  (u),  is  chosen  for  the  average.  The  cancellation  of  this  term  and  the 
oscillatory  term,  ),  when  m  =  n  also  depends  on  overlap  of  the 

receiver  illumination  functions  dj^  and  d^;  hence  the  desirability  of  this  con¬ 
dition. 

(C)  The  term,  introduces  additional  delays,  to  the  target- induced 

delays,  However,  these  delays  are  expected  to  be  small  compared  with 

significant  target- induced  delays  at  large  ranges  for  the  following  reasons. 
Reduction  of  data  taken  by  Mackenzie  [11],  and  by  Nuttal  and  Cron  [12]  indi¬ 
cates  that  0.  2>  c^  >  0.  06  and  0.  2  <  <  1  for  angles  of  incidence  in 

the  range  60  <  Y  s  90°.  Also,  penetration  depths  of  less  than  2  m  are  expected 
for  these  angles  of  incidence  [  12] .  Assuming  that  the  reflection  coefficient 

model  of  Equation  (A  14)  holds,  and  that  no  more  than  five  internal  reflections 

r  /O  2 1 

within  the  top  sediment  layer  are  significant,  p^  <  (5)  (-—1(2)  ~  1  msec  might 
be  expected.  However,  the  discussion  of  Section  II  indicates  that  significant 
target  delays  measurable  by  STARLITE  techniques  are  greater  than  10  msec, 
roughly. 

2  2  2 

(C)  The  condition,  Q  «  A  «  u^,  ensures  that  u  -  n  /4  w  .  Thus,  the 

essential  u  dependence  of  the  diagonal  terms  in  Equation  (64)  reduces  to 

2  i 

functions  of  the  form  u  exp  (-2wp  ).  Recall  tha^  |3  =  0,  by  definition;  hence, 

2  ^  °  i  - 1 

the  first  term  in  the  sum  is  ~  u  while  terms  for  n  s:  1  have  maxima  at  u  =  (p„  ) 

For  0.2  sc  i3,  s  1,  these  maxima  occur  at  frequencies  to  <  7.  5  x  10  sec 
Thus,  at  carrier  frequencies  of  f^  >  1  KHz,  the  spectrum,  I(u,  n).  is  certainly 
not  flat.  Therefore,  it  is  necessary  to  rely  on  the  condition  A  «  to^  to  pre¬ 
serve  a  semblance  of  flatness  for  I(u,  n)  in  the  integration  interval  [w^  -  A/ 2, 
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(C)  With  regard  to  the  off-diagonal  terms,  m  7^  n,  in  Equation  (64),  it  is 
reasonable  to  expect  that  if  exp  [l/2fi(p^^  -  is  significantly  different 

from  unity,  then  exp  much  smaller  than  a  comparable 

diagonal  term  exp  [-2w|J.  ^],  and  therefore,  such  off-diagonal  terms  are 
negligible.  Off-diagonal  terms  that  are  not  negligible  in  comparison  with 
like  diagonal  terms  require  that  )  «  1  in  order  that  I(u,  fi)  not 

be  too  oscillatory  over  the  integration  band. 


Rough  Surfaces 


(C)  Next  we  consider  the  effect  of  surface  roughness  by  assuming  that 

2 

w  no  »  1,  but  that  the  a  are  sufficiently  small  to  neglect  the  K  .  Then, 

®  2  C  2  2 

Vm  ^  ^  exponential  integral  in  takes  on  its  asymtotic 

value  and  dominates  the  logarithmic  term.  If  (L/D)^  «  1  by  virtue  of 

short  correlation  lengths  L, 

2  2 


I(w,  Q )  const 


If  (I  /D)^  E.  »  1, 


I(u>  Q)  const  •  e 


"n  CT. 

2  C 

(j  e 
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,  a 
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n,  m=o 


X.  \  e 
n  m 
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01 

1  2/  m 

(66) 


1  2  2^2 

2  h  a 
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n,  m=o 


u 


.Q' 


+  3  Iw+f' 


m 


\  \ 
n  m 


(67) 


2  2  2 

In  either  case,  r|  o^  Q  «  1  is  required  if  a  strong  amplitude  dependence  on  Q 

is  to  be  avoided.  Measurements  of  a.  in  the  Hatteras  Abyssal  plain  indicate 

s 

that  <  0.  5  m  in  this  region  [13].  At  incident  anglec  f  a  60°  and  s  0.  5  m, 

2  2^2  ^ 

Q  Oq  t]  s  0. 1  for  Q  s  900.  This  is  well  within  the  range  of  frequency  shifts 

for  which  the  STAR  LITE  method  is  applicable.  It  is  also  to  be  noted  from 

Equations  (66)  and  (67)  that  the  u  dependence  of  I(u,  n)  is  strongly  affected  by 

the  surface  roughness  and  its  correlation  lengths  relative  to  the  dimensions 
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(C)  of  the  illuminated  surface.  In  fact,  rough  surfaces  with  long  correlation 
lengths,  Equation  (67),  more  closely  approximate  the  desired  flat  frequency 
spectrum  of  I(u,  n)  than  do  smooth  surfaces.  Equation  (64). 


Irregular  Sediment  Layers 


(U)  Lastly,  we  consider  briefly  the  effect  of  irregularities  in  the  sediment 
layer  thicknesses  on  I(u,  Q).  Defining 


K  =  K  -  iK  '  =  L  (/3  ^13/ V  ^  -j  L  (|3  '  +  i3"')cj  I  (68) 

"  "  "  q,p=o  ^  '^q'^p'qp  ''q,p=o  ‘I  P  P  ‘IP 


the  exponent  of  interest  is  of  the  form 


X  =  - 


2 

UQ 


w-8-1  K  - 


h  2 


U+f  K 
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r 


n  -  ^  r^2 


1 

2 


r  v 
K  +  K 
n  m 


-  wo 


K  ^  +  K  ^ 
n  m 


(69) 


Again,  consider  the  case  of  incident  angles  near  and  above  critical,  and 

assume  that  essentially  only  the  top  layer  is  penetrated,  so  |3^  =  IS for  all  n. 

Near  critical  incidence,  ^  IS.^,  so  K  ^=“0.  However,  abovecriticalin- 
j  1  1  n  2  r 

cidence,  ^  n  creates  terms  of  the  form  exp(+  w  |K  | 

in  I(w,  0)  which  tend  to  compensate  for  the  negative  slope  of  terms  of  the  form 

exp(-W|J.  V  The  terms  l/2wn()^^  -  K^^)  and  1/2  j  (w^  +  a^/4)(K^^  -  K^^)  in 

X  drop  out  in  the  diagonal  elements  of  the  double  sum. 


2  r 

(C)  To  estimate  the  effect  of  the  term,  1/4  n  |  |,  in  X,  consider  inciident 

angles  near  grazing  and  assume  that  c^^^^  =«  1,  ^  0.06/c^,  and  that  no 

more  than  five  internal  reflections  are  significant  as  before.  Then,  if  n  ^ 
600  sec’ 1/4  jKg^l  5  0. 1  if  <  0.3  m.  This  value  of  is  compar¬ 
able  with  the  value  of  quoted  for  the  Hatteras  Abyssal  Plain.  It  therefore 
appears  that  if  the  surface  is  not  too  rough  for  STARLITE,  neither  will  the 
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(C)  sediment  layers  be  too  irregular.  The  remaining  term,  -1/2  3wn(K^^+ 
in  X  may  be  troublesome,  however,  since  under  the  above  conditions  of  near 
grazing  incidence,  <4  x  10  w,  which  is  not  necessarily  small  for 

carrier  frequencies  of  f^  >  1  KHz. 

(C)  From  this  brief  and  somewhat  superficial  discussion  of  the  properties  of 
I(u,  Q),  it  seems  reasonable  to  conclude  that  there  is  some  justification  for 
believing  that  STARLITE  techniques  are  applicable  to  signal  classification  in 
the  bottom  bounce  mode  of  operation.  The  conditions  for  operation  at  long 
range  appear  to  be  most  favorable,  since  the  effects  of  surface  roughness 
and  penetration  of  the  sediment  layers  by  the  acoustic  wave  are  reduced  at 
angles  approaching  grazing  incidence.  However,  the  degree  of  confidence  in 
success  that  one  would  like  before  initiating  costly  sea  trials  for  verification 
of  the  classification  method  has  not  been  achieved  in  this  study.  More  exten¬ 
sive  analysis  of  the  properties  of  Equation  (55),  as  a  function  of 

geometry,  signal  spectrum,  and  parameters  of  the  bottom  transfer  functions 
through  computer  simulation  is  expected  to  provide  a  higher  confidence  esti¬ 
mate  of  the  probability  of  success  of  bottom  bounce  STARLITE.  In  addition, 
the  effects  of  target  and  transmitter-receiver  motion,  ocean  surface  back- 
scatter,  volume  and  surface  reverberation,  ambient  noise,  and  scattering 
off  corrugated  surfaces  should  also  be  studied  before  final  judgment  is  passed 
on  the  theoretical  practicality  of  bottom  bounce  STARLITE. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


(C)  To  establish  direction  to  this  study  of  souar  signal  characteristics  in  the 
bottom  bounce  mode,  it  was  decided  to  concentrate  on  the  effects  of  bottom 
reflections  as  related  to  target  classification  on  the  basis  of  shape  and  aspect 
for  targets  characterized  as  arrays  of  point  scatterers.  A  review  of  the  gen¬ 
eral-theoretical  requirements  for  target  shape  and  aspect  classification 
through  spectral  and  time  crost.- correlations  of  signals  received  at  two  dif¬ 
ferent  space  points  has  indicated  that:* 

1)  Spectral  processing  is  more  promising  than  time  processing 
for  shape  and  aspect  classification  from  a  practical  point  of 
view,  since  considerably  larger  receiver  separations  are  re¬ 
quired  for  time  correlations  (of  the  order  of  the  ratio  of  the 

signal  carrier  frequency  to  the  bandwidth,  ff,/B  ,  times  larger). 

u  w 

We  identify  the  special  case  of  time  cross-correlation  of  the 
envelopes  of  uncompressed  linear  FM  signals  with  spectral 
processing  in  this  observation. 

2)  A  band-limited  spectral  cross-correlation  with  symmetric 
frequency  shifts  of  the  form 

,  A 

R^2(^)  =  Ref  ^  W^(u)  -  +^)  du 

"  J 

■'•The  authors  regret  that  they  were  unable  to  obtain  a  copy  of  the  report  "An 
Experimental  Study  of  STAFtLITE  using  a  Scaled  AN/SQQ-23  (PAIR)  Config¬ 
uration,"  DRL-TM-68-10  (AD  390  245L),  in  time  to  influence  this  theoret¬ 
ical  study. 
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should  be  used  for  spectral  processing  to  confine  the  ob¬ 
servation  of  the  target  transfer  function  to  the  interior  of  the 
spectral  window  created  by  the  transmitted  signal  and  ne¬ 
gate  the  errors  in  estimating  target  aspect  introduced  by 
attenuation  in  bottom  reflections. 

3)  Geometric  and  signal  parameter  constraints  for  spectral 
processing  are  somewhat  different,  and  more  restrictive 
in  the  case  of  the  correlation  and  ambiguity  conditions,  but 
less  restrictive  in  the  case  of  the  resolution  and  transient 
conditions,  than  indicated  in  the  original  STARLITE  analysis 
[2],  [23]  of  spectral  processing. 

4)  Large  f^  and  are  desirable  for  satisfying  spectral  pro¬ 
cessing  constraints  with  reasonable  receiver  separations, 

whereas  low  and  small  B  are  desirable  to  minimize  prop- 
u 

agation  loss  and  distortion  of  the  spectral  window  in  bottom 
reflections.  Therefore,  optimum  f^  and  depend  upon  tar¬ 
get  range  and  bott.  m  conditions. 

5)  An  example  for  f^  =  6  KHz  and  B  ^  =  600  Hz  indicated  that  re- 

u  w 

ceiver  separations,  D,  of  about  175  yards  permitted  class¬ 
ification  over  line  target  aspect  angles  of  15°  ^  cd  ■s  60°  and 
120°  ^  Cv  ^165°  (bear  sp'^ct  at  0°  and  180°,  bow  or  stern  as¬ 
pect  at  90°  in  this  nocdiion)  for  a  target  of  length  =  90  m 
at  a  bearing  9  =  90°  and  ranges  of  20  to  40  thousand  yards. 
Howev^er,  since  D  is  inversely  proportional  to  L^,  f^,  and 
sin  9,  reduction  of  any  of  these  parameters  requires  a  cor¬ 
responding  increase  in  D  to  maintain  classification  over  the 
same  interval  of  ranges  and  aspect  angles.  It  therefore 
appears  that  some  flexibility  in  selecting  receiver  separations 
as  a  function  of  range  and  bearing  is  necessary. 
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(C)  V/ith  regard  to  the  feasibility  of  STARLITE  processing  in  the  bottom  bounce 
lode,  our  conclusions  are  somewhat  more  tentative.  A  physical  approach  to 
the  modeling  of  sonar  wave  reflections  off  an  irregular"  impedance  boundary 
[7]  has  been  taken  ir  the  analysis  in  order  that  the  frequency  dependance  of 
the  bottom  transfer  function  can  be  related  to  measurable  physical  properties 
of  the  bottom  sediments.  The  analysis  requires  a  considerable  degree  of 
approximation  and  a  number  of  assumptions  to  obtain  analytically  tractable  re¬ 
sults.  However,  comparisons  of  similar  analytic  developments  with  con¬ 
trolled  experiments  indicate  that  the  assumptions  and  approximations  are  per¬ 
haps  not  too  restrictive  as  applied  to  this  type  of  problem  [?]: 

1)  STARLITE  signal  processing  concepts  indicate  that  the  ideal 
bottom  transfer  function  would  be  independent  of  frequency 
within  the  interval  chosen  for  spectral  cross-correlation. 

However,  it  has  been  shown  that  the  spectrr'  properties  of 
the  bottom  transfer  function  can  vary  widely  (low-pass,  high- 
pass,  bandpass,  or  notched-filter  shape)  depending  upon  L,  B.  , 
angle  of  incidence,  the  degree  of  surface  roughness,  irregularities 
in  the  sediment  layer  thicknesses,  and  the  amount  of  attenuation 
associated  with  the  sediment  layers. 

2)  In  spite  of  the  complexity  of  the  bottom  transfer  function,  a 
cursory  examination  of  its  frequency  dependence  for  para¬ 
meter  values  typical  of  the  abyssal  plains  leads  to  the  con¬ 
clusion  that  long-range  bottom  bounce  STARLITE  could  be 
practical,  since  the  effects  of  surface  roughness  and  pene¬ 
tration  of  the  boundary  by  the  acoustic  wave  are  minimized  as 
the  incident  angle  approaches  grazing  incidence. 

(U)  Because  of  the  large  number  of  parameters  which  characterize  the  class¬ 
ification  problem  in  the  bottom  bounce  mode,  and  the  uncertainity  in  their 
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values  and  distributions,  a  parametric  study  through  computer  simulation  of 
the  classification  scheme  is  a  recommended  step  before  the  expense  of  veri¬ 
fying  bottom  bounce  STARLITE  classification  techniques  through  special  sea 
trials  could  be  justified.  Such  a  simulation  would  not  require  the  many  assump 
tions  and  approximations  required  in  the  present  study  and  would: 

1.)  Establish  theoretically  whether  or  not  bottom  bounce  STARLITE 
classification  concepts  are  feasible. 

2)  Provide  valuable  guides  to  possible  future  experiments  at  sea 
by  determining  optimal  signal  carrier  frequencies  and  band- 
widths,  and  parameter  constraints  that  are  realistic  to  the 
bottom  bounce  mode  of  operation. 
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APPENDIX  A 

BOTTOM  REFLECTION  COEFFICIENT  MODEL 


The  basic  model  for  the  ocean  bottom  is  a  regular  stratified  medium  as 
shown  in  Figure  Al,  All  layered  media  are  assumed  to  be  absorbing,  and 
for  the  moment,  separated  by  plane  boundaries.  The  bottom  is  assumed  to 
consist  of  N  layers  of  finite  thickness,  -t^,  1  ^  n  ^  N,  overlying  an  absorbing, 
semi-infinite  supporting  half-space  (medium  N  +  1).  The  ocean  is  defined  to 

be  the  o^"  medium  and  is  assumed  to  be  a  semi-infinite,  nonabsorbing  fluid. 

th  *  th 

We  define  the  interface  to  be  the  boundary  between  the  n  “  and  (n  +  1) 

media. 

From  wave  impedance  concepts  for  plane  waves  with  harmonic  time 
dependence  the  reflection  coefficient  at  the  n^^  interface  for  obliquely 

incident  waves  is  given  by  [14]; 


^  4.  1  e 

-_I1 _ n.t-l 


n 


1  +  X.  A  .16 
n  n  +  1 


— jf 


(Al) 


n  +  1  n  +  1 


til 

In  Equation  (Al),  A  is  interpretable  as  the  reflection  coefficient  at  the  n 

^  Q+ 

interface  for  waves  impinging  on  the  (n  +  1)"  to  (N  +  1)  media  when  the  n 

A 

medium  is  considered  to  be  semi- infinite,  X  is  the  reflection  coefficient  at 

+h  th  ^st 

the  n  interface  when  both  the  n  and  (n  +  1)  media  are  considered  to  be 

A  1 

semi-infinite,  and  jS  is  one-half  the  normal  component  of  the  complex 

^  th 

phase  velocity  of  the  wave  propagating  in  the  n  layer. 


To  relate  X^^  and  to  acoustical  properties  of  the  media,  consider,  for 
example,  plane  wave  propagation  in  an  absorbing  fluid  for  which  the  Naviar- 
Stokes  and  continuity  equations  take  the  following  form: 
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PN+1,  '^N+l,  ‘‘N+I 


Figure  Al.  Geometry  for  a  Regular  Stratified  Bottom 
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3- 

"Sf  '  “ 

t  n 


+  p  div  u 
n  n 


0 


(A2) 


In  these  equations,  is  the  pressure,  the  vector  particle  velocity,  the 
fluid  density,  r^  the  normalized  acoustic  flow  resistivity,  and  c^  the  phase 
velocity  in  the  limit,  r^  -  0.  The  intrinsic  impedance,  z^,  of  a  medium 
characterized  by  Equation  (A 2)  is 

\  ’  "n  ‘^nV  ‘  j 

and  the  complex  wave  number  for  harmonic  plane  waves  is 


k  =  j  r  /u 

n  c  V  n 

n 


(A4) 


For  a  plane  wave  incident  from  the  0  medium  at  a  (real)  velocity,  c^,  and 
angle,  Snell' s  law  yields 


k  sin  Y  =  k  sin  'f 
n  n  o  o 


(AS) 


where  k  =  w/c  .  From  Equations  (A4)  and  (AS),  the  cosine  of  the  (complex) 
o  o  th 

angle  of  incidence,  at  the  n''"  interface  is 


(A6) 
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Since  the  complex  phase  velocity  in  the  n 


th 


medium  is  from  Equations 


(A4)  and  (A 6)  and  the  definition  of 

^  2  cos  Y  /(I  -  j 


2 

n 


.  2  V 

SUl  I, 


(A  7) 


Finally,  in  terms  of  the  normal  wave  impedance,  z^,  the  complex  reflection 
coefficient,  is  given  by  [l4] 


A 

\ 

n 


"n  +  1 


^n  +  1 


+  z 

n 


where 
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Jl. 


cos  Y 

n 


2  p  (1  -  j  r  /u) 
n  _ y  _  a  ■ 


^n 


(A8) 


(A9) 


For  future  reference,  we  note  that  physically,  r  >  0.  Thus,  by  Equation 

A 

(A  7),  /3^  is  of  form 


0  =  0  \  0  ^>0,  0  ^>0 


(A  10) 


for  positive  frequency,  u;  >  0,  0  and  any  finite  ratio  of  velocities. 


The  recursion  relation.  Equation  (Al),  terminates  on  the  condition, 

A 

principle,  can  be  used  to  generate  a  closed-form  solution 
for  the  desired  reflection  coefficient  at  the  bottom,  A^.  However,  in  practice 
we  shall  have  to  be  content  with  a  series  expansion  of  A  carried  to  some 
acceptable  order  of  accuracy  in  powers  of  or  exp  (-  u0^  -t^).  To  deter¬ 
mine  the  nature  of  these  expansions,  consider  briefly  some  properties  of 


46 


UNCLASSIFIED 


UNCLASSIFIED 


bottom  sediments  and  the  bottom  reflection  coefficient  as  determined  from 
bottom  reflectivity  experiments  [l5,  16],  Measurements  of  the  acoustic 
properties  of  bottom  sediments  indicate  that  the  sediments  near  the  water 
interface  have  nearly  the  same  speed  of  sound  as  water  and  that  the  velocity 
and  density  tend  to  increase  with  depth.  Adjacent  sediment  layers  seem  to 
have  nearly  the  same  intrinsic  acoustic  impedance,  so  that  for  near  normal 
incidence,  <<  1.  A  series  expansion  of  Equation  (Al)  is  of  the  form 


n 


i+h  ,  .e'j'^^n  +  I'^n  +  1 

n  n  +  1 


E 

m=0 


<-^n  ^n  +  P  ^  +  l  n  +  1 


(All) 


which  indicates  that  to  second-order  in  X  A  is  linear  in  X  .  That  is, 

n’  n  n 


A  -  X  +  A  e"  j“^n  +  l^n  +  1 
n  n  n  +  1 


^e'^'  X  e-j^"  S  ^n  +  q^n  +  q 

k=o  q=o 


(A  12) 


and,  in  particular, 


N 


n 


A  (2)  .  s  ;  e-  ■*“  S 
°  n=0  ”  q=0  ‘I  1 


(A  13) 


A 

where  jS  =  0.  Bottom  reflectivity  measurements  in  [l6j  indicate  that 
A  <  0.  4  for  indident  angles  of  <  60^.  Since  the  acoustic  wave  penetration 
of  the  boundary  was  apparently  quite  deep,  and  encompassed  several  acousti- 

A 

cally  distinct  sediment  layers,  one  expects  that  the  X^  for  these  measure¬ 
ments  were  somewhat  smaller  than  A^  and  that  the  approximation  of  Equa¬ 
tion  (A  13)  would  hold  over  this  range  of 
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Measurements  of  the  reflection  coefficient  [16]  over  the  range  60®  <  'i'^<  90' 


indicated  that  the  penetration  of  the  bottom  by  the  acoustic  wave  was  quite 
shallow  (<  2  m).  If  we  assume  that  essentially  only  one  sediment  layer  is 
penetrated  at  angles  near  and  above  critical  incidence,  from  Equation  (Al), 


with  Aj  = 


^o  =  ^  ® 

n=o 


(A  14) 


The  series  in  Equation  (A14)  can  be  terminated  after  N  terms  when 
is  sufficiently  large  and  the  equation  written  in  the  same  form  as  Equation  (A13). 


1.  e. , 


n  ~ 

A,  L  ^  I 

N  ^  q=o 
\  £  A  e 


o  n 

n=o 


(A15) 


where 


(A16) 


0;  0  =  |3.,  q  =  1,  . . .  ,  n 
Hi  ^ 


For  the  purposes  of  the  discussions  in  this  report,  then,  it  will  be  assumed 
that  the  bottom  reflection  coefficient.  A,  is  adequately  described  by  a  function 
of  the  form 


n  q=  O  "^q 


(A  17) 


where  (3^  =0,  or  and  ,"1^  =  or  3^,  whichever  is  appropriate  to  the 


angle  of  incidence, 
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Since  a  is  an  important  contributer  to  the  transfer  function  representing 
bottom  reflections,  the  frequency  dependence  and  randomness  of  the  param¬ 
eters  describing  A  are  of  interest.  With  regard  to  the  frequency  dependence, 
measurements  of  attenuation  in  bottom  sediments  indicate  that  it  is  a  linear 
function  of  frequency  [  12,  17,  18,  19].  For  fluid-like  sediments  describable 

A 

by  Equation  (A2),  this  implies  that  the  j3  are  independent  of  frequency,  and 
hence  that  =  constant.  Therefore,  by  Equations  (A8)  and  (A9),  the 
are  also  independent  of  frequency.  With  regard  to  randomness  in  the  reflec¬ 
tion  coefficient  over  the  area  illuminated  by  the  transmitter  and  receiver, 
inhomogenities  in  the  sediments  can  create  variations  in  the  density  and  phase 
velocity  of  the  media,  where  uneven  deposition  of  sediments  gives  rise  to 
variations  in  the  layer  thicknesses  and  mixing  of  sediments  at  the  boundaries, 
thus  making  it  difficult  to  distinguish  a  precise  boundary.  Therefore,  it  is 
expected  that  the  quantities  j3^,  and  are  best  described  by  random 
functions  of  the  spatial  coordinates  of  the  bottom,  but  not  time.  However, 
for  simplicity,  v/e  will  ignore  any  possible  randomness  ii^  thereby 

assuming  that  the  essential  nature  of  the  randomness  in  A  is  contained  in  the 
layer  thicknesses,  In  summary,  the  model  chosen  for  the  bottom  reflec¬ 
tion  coefficient  is  that  of  Equation  (A17),  where  X^  and  are  complex  con¬ 
stants  and  t  is  a  real,  random  function  of  the  bottom  coordinates. 
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APPENDIX  B 

EVALTTATION  OF  THE  INTEGRAL,  I(u,  fi) 


In  this  appendix,  we  consider  a  solution  of  the  integral 


Ku.  n)  =  I 

4  /  y  y 


-jTif|u> -f)c(x)  -  + 

d0d(x')<Aj^.u a"'  ;.'u  +  f><e  : 


/  (Bl) 


under  the  assumptions  that: 

1)  The  bottom  apertures  are  gaussian  illumination  functions: 
d{x)  =  d(x,  y)  =  exp 


,2  2  n 

1 1 

"2 

Dx^  Dyl 

(B2) 


2)  The  reflection  c  efficients  have  the  form  stipulated  in 
Appendix  A: 


N 


-ju  ^  t  (y) 
■’  q=o  *^0  q 


A(;u>)  =  S  X^e  ■'  q 


n=o 


(B3) 


3)  The  sediment  layer  thicknesses,  I  (y),  and  surface  roughness, 
C(y),  are  gaussian  random  variables  with  homogeneous  correla¬ 
tion  coefficients.  More  explicitly,  for 


A  A 


X  -  y '  =  (x  -  X  0  i  +  (y  -  y ')  3  =  ui  +  vj 


(B4) 
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it  is  assumed  that 


<C(x)>  =  0 

<C(x)  C(x')  >  =  v) 

=  <Q^(x)  > 

p^(0,  0)  =  1 


(B5,) 


and 


“qp  Pqp<“'''> 

p  =  P  .  p  (0,  0)  =  1 

qp  pq  qp 


(B6) 


Osa  ^  -  <(i  -  I  )^> 

qp  pq  q  q 

2 

In  the  notation  of  Equations  (B5)  and  (B6),  a  is  a  variance  and  p  is  a 

2  2  . 

normalized  spatial  correlation  coefficient.  The  variance,  a  ^  a  ,  is 
specified,  since  different  sediment  layer  thicknesses  at  the  same  place  do 
not,  in  general,  correlate  perfectly. 

Denoting  C(x')  and  by  C'  and  respectively,  l(u,  u)  takes  the 

form 


Kcj,  n)  = 


2  Q" 

-T 


N  N 
S  2  X  \ 
n=o  m=o  ^  XX 


d(x)d(x')  <e 


-1 


M 


n 


<e 


•  I 

-jqju  -  ^ 

n 

E  IS  I 

q=0  “I 


c  - 


\> 


m 

L  IS  "  I  'i> 

p=o  P  P 


(H7) 
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It  is  well  known  that  the  joint  characteristic  function,  Y^(v).  of  a  vector 
gaussian  variable,  z,  z  =  [z^,  . . . . ,  z^],  is  [20]: 


.  t 

Y  (v)  =  <e"^-  ->  =  e 
z 


.  t-  ^  t  ^ 
-jv  z  -  2v  K  y 


(B8) 


where  v  is  an  n- dimensional  column  vector,  t  denotes  the  transpose,  z 
is  the  mean  of  z,  and  K  is  the  correlation  matrix 

K  =  <(z  -  z)  (z  -  z)*> 

-  -  -  - 


(B9) 


Since  every  expectation  in  I((j,  fi)  is  of  the  form  of  Y,  we  can  write 


I((J,  Q)  = 


^ "x?  f  <B10) 


,2\  N 


The  characteristic  functions,  and  Y,  ,  are  determined  simply  from 


— nm 


Equation  (B8)  and  the  definitions  of  u>  K  ,  u  >  K  .  once  ^  and  L 
^  ^  g  — nm  '"_nm  —  — nm 

are  defined. 


If  we  define  the  two-dimensional  vectr.-r  of  surface  roughness  parameters  to 
be 


C  =  K-S'] 

from  Equation  (B7),  we  recognize 


t 


10-9 

2 


From  the  assumed  properties  of  Q  in  Equation  (B5) 

1  PC 


Sc  =  "c 


PC 


C  =  0 


hence 
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where  the  K  are  the  (2  x  2)  matrices 

~qp 


K  =  a 
-v<qp  qp 


qp 


^qp 


=  K 


pq 


and  is  the  (n  -  m)  x  (n  -  m)  matrix  with  elements 


m  +  i,  m  +  j 
Then,  noting  that 

u  ^  u  =  a  ^ 

-q  ,^qp~p  qp 


Isi,  j  5n-m 


L  -  fl^|3  ^  --^1(0  *  +  13  |3  )p 

](  2 1  I  T  p  p  '^q'^pj'^qp 


(0 

2  .,1 
(3  3  1 

2 

q  p  f 

) 

J 

=  e 

inm 


CO 


n 

L  |3  I  +  j 
\2  n 


0 


m 


^  q=o  ^ 


w-§)l  130^ 

2)qp=o  'iP^P  2 


(B12) 


o  ^21  Q 
2  O  nm  2 

^  I  L  Re 

q,p-o  'iP 


q  p 


w+f  ?  0  “'a  ^ 

2  L  q  ^p  qp 

q,p=o 


Pqp  <“• 


where  =  min(n,  m). 


Consider  next  the  typical  spatial  integral  in  Equation  {B7)  by  picking  out  only 
those  terms  in  Equations  (Bll)  and  (B12)  which  are  functions  of  y  and  y'. 
The  typical  integral  is 
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■^nm"  OT  dxdydx'dy^e 


1 1  x^^  ^ 

Dx^  Dy^ 


(B13) 


Define 


'  2  Q 

(j 


Q- 


2  2 


Q 


nm 


Ti  Or  (u.v)+  L  o  Re 
'•  q.p=o 


p„„(u,  v) 

q  p 

^qp 

2  2  2  nm  2 

t  =Tia-+L  ^  Re 

"""  ^  q.p=o. 


q  p 


(B14) 


and  a  composite  c'^rrelation  coefficient 


o'  (u,  v)  =  7  J 
'^nm'  t 


nm 


Q 


.y  9  nm  9 

T1  o  p  (u,  v)  +  L  C  Re 
^  ^  q,  p=o 


^q'^plPnm 


(B15) 


Then,  substituting  x=x'  +  u,  y=y^+v  and  performing  the  x*  and  y*  inte- 

2\ 


gration,  reduces  to 


If  u'  v" 


rr  ^'Dx^  Dy^ 

Jr.rv,  “  TjDxDy  Jj  dudv  e  ^ 


0,2. 

4 


nm 


(B16) 


Equation  (B16)  cannot  be  integrated  exactly  for  the  usual  forms  of  How¬ 

ever,  certain  reasonable  functions  permit  approximate  solutions  in  terms  of 

common  tabulated  functions.  With  this  goal  in  mind,  it  is  assumed  that  p 

°  ^nm 

is  of  the  form 


L^^nm 


Ly 


nm 


2  a  uv 
nm 

Lx  Ly 
nm  nm 


(B17) 


where  Lx^^^  and  Ly^^  are  correlation  lengths  appropriate  to  the  principal 

axes  of  d(v )  and  0  s  a  <1,  with  a  =0  when  the  principal  axes  of  S'  „ 
^  nm  nm  '^nm 

coincide  with  those  of  d(^).  Expanding  the  exponent 


55 


UNCLASSIFIED 


UNCLASSIRED 


L 

k=o 


k! 


'run 


(B18) 


n 


* 

- 

I 

J  _  n^J 

t 

,,2  (I2l  2 
T  rnm 

k  i  '  ; 

4  1 

nm  '^nm  <» 

/  J 

^  W-  i 

J  can  be  written  as  the  series 
nni 


run 


00 

L  J 
k=o 


nmk 


(B19) 


where  J  ,  is  of  the  form 
nmk 


TtDxDy 


“  4 


nm 


runk 


k! 


!  ^5 


1  t^  -1 

2  “  <,^0  - 


t 


with  z  =  [u,  V  ]  and 


K 


1 


2Dx 


Lx^  / 
run  ' 


ka 


nm 


Lx  Ly 
nm  nm 


ka 


nm 


Lx  Ly 
nm  •'nm 


1  +  k 


2Dy^  Ly^ 

run 


(B2a) 


(B21) 


Equation  (B20)  is  of  the  form  of  a  bivariate  gaussian  integral  with  the  result 
that 


2tt  DxDy 


4  I  run 


nmk 


k!  |K  "^1 
'~o  ' 


(B22) 


For  bottom  reflections  far  from  the  transmitter  and  receiver,  it  is  reason- 

2  2  2  2 

able  to  assume  that  Dx  /Lx  »  1  and  Dy  /Ly  »  1.  Hence, 

/  nm  •'  '  run 


l5o''  I  “< 


k  =  0 


2DxDy  ' 

_ ./t: 

Lx  Ly  V 
nm  nm  ’ 

56 


a  k  a  1 

run 
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